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ff^ ' Abstract 

Chains with long-range interactions are considered. The interactions are defined such 
f^ '_ that each nth particle interacts only with chain particles with the numbers niba(m), where 






X 



m = 1, 2, 3, ... and a{m) is an integer-valued function. Exponential type functions a{m) 
6"*, where h = 2,3,.., are discussed. The correspondent pseudodifferential equations of 
chain oscillations are obtained. Dispersion laws of the suggested chains are described by 
the Weierstrass and Weierstrass-Mandelbrot functions. 
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Long-range interaction (LRI) has been the subject of investigations for a long time. An infinite 
Ising chain with LRI was considered by Dyson pQ. The d-dimensional Heisenberg model with 
long-range interaction is described in Refs. [21 E], and its quantum generalization can be 



?H ' found in lU [5]. Solitons in a chain with the long-range Lennard- Jones- type interaction were 

considered in [7]. Kinks in the Frenkel-Kontorova chain model with long-range interparticle 
interactions were studied in [S]. The properties of time periodic spatially localized solutions 
(breathers) on discrete chains in the presence of algebraically decaying LRI were described in 
[9]. Energy and decay properties of discrete breathers in chains with LRI have also been studied 
in the framework of the Klein-Gordon [6], and discrete nonlinear Schrodinger equations [10]. 
A main property of the chain dynamics with power-like long-range interactions [9l [11] is that 
the solutions of chain equations have power-like tails. The power-law LRI is also considered in 
Refs. [121 [121 HH [I5]- In Ref. [21], we formulate the consistent definition of continuous limit 
for the chains with long-range interactions (LRI). In the continuous limit, the chain equations 
with power-like LRI give the medium equations with fractional derivatives. 



Usually we assume for LRI that each chain particle acts on all chain particles. There are 
systems where this assumption cannot be used. In general, the chain cannot be considered as 
a straight line. For example, the linear polymers can be represented as some compact objects. 
It is well known that "tertiary structure" of proteins refers to the overall folding of the entire 
polypeptide chain into a specific 3D shape [161 [13 [IE]- The tertiary structure of enzymes is 
often compact, globular shaped [HI [I?]. In this case, we can consider that the chain particle is 
interacted with particles of a ball with radius R. Then only some subset A^ of chain particles 
act on nth particle. We suppose that nth particle is interacted only with fcth particles with 
k = n± a{m), where a{m) G N and m = 1,2, 3, .... We can consider fractal compactified linear 
polymers (chains), such that these "compact objects" satisfy the power law N{R) ~ R^, where 
2 < D < 3 and N{R) is a number of chain particles inside the sphere with radius R. As an 
example of such case, a{m) will be described by exponential type functions aim) = 6"*, where 
6 > 1 and 6 G N. In this case, the LRI will be called fractal interaction. 

The goal of this paper is to study a connection between the dynamics of chain with fractal 
long-range interactions (FLRI) and the continuous medium equations with fractal dispersion 
law. Here, we consider the chain of coupled linear oscillators with FLRI. We make the transfor- 
mation to the continuous field and derive the continuous equation which describes the dynamics 
of the oscillatory medium. We show how the oscillations of chains with FLRI are described by 
the fractal dispersion law. This law is represented by the Weierstrass functions whose graphs 
have non- integer box-counting dimension, i.e., these graphs are fractals. Fractals are good mod- 
els of phenomena and objects in various areas of science [20]. Note that fractals in quantum 
theory have recently been considered in [21]. In this paper, we prove that the chains with 
long-range interaction can demonstrate fractal properties described by fractal functions. 

2 Chain equation 

One of the oldest fractal functions is Weierstrass function [22] : 

OD 

ly (x) = ^ 6" cos(a"7rx) (1) 

n=0 

introduced as an example of everywhere continuous nowhere different iable function by Karl 
Weierstrass around 1872. Maximum range of parameters for which the above sum has fractal 
properties was found by Godfrey Harold Hardy [23] in 1916, who showed that 

< 6 < 1, ab> I. 



The box-counting dimension of the graph of the Weierstrass function W{x) is 

(2) 
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Functions whose graphs have non-integer box-counting dimension are called fractal functions. 
Consider a one-dimensional system of interacting oscillators that are described by the equa- 
tions of motion, 

— ^M„(t) = C^ 22 ~~p2~ \^n+a{ni){t) - '2Un{t) + M„_a(m)(t)], (3) 



m=—oo 



where Un{t) are displacements from the equilibrium, and h is the distance between the oscilla- 
tors. Here a(r7i) and 6(m) are some functions of integer number m. The case b{m) = J{\n — m\) 
has been considered in [24|. The right-hand side of the equation describes an interaction of the 
oscillators in the system. 

We illustrate this chain equation with well-known example [26]. In the case of nearest- 
neighbor interaction, we have 

a{m) = 1, 6(m) = 6mo- 

Then equation (^ gives 

^2 2 

— u„(t) = -^ K+i(t) - 2n„(t) + u„_i(t)]. (4) 

We can define a smooth function u{x,t) such that 

u{nh,t) = Unit). (5) 

Then equation (JH) has the form 

d^u{x, t) = -— [u{x + h,t) — 2u{x, t) + u{x — /i, t)]. (6) 

This is the differential-difference equation. Using the relation 

exp{—ihd^) u{x, t) = u{x + h,t), 
we can rewrite equation ([6]) in the form 

dMx, t) + -^ sin^ (- Y^,) u{x, t) = 0. (7) 

This is the pseudo-differential equation. The properties of this equation have been considered 
in [2S]- For h —^ 0, we obtain the wave equation 

d'}u{x,t) -c^dlu{x,t) =0. 



3 Fractal long-range interaction 

If a{m) in equation ([1]) is not a constant function, then we have the long-range interaction of the 
chain particles. Note that the function a{m) should be integer-valued. For example, a{m) = m 
and a{m) = 2™. The set A^ = {n±a{m) : m G N} describes the numbers of particles that act 
on the nth particle. 

(a) If a(m) = m, where m G N, then An is a set of all integer numbers Z for all n, i.e.. 
An = Z. In this case, the nth particle interacts with all chain particles. 

(b) If a{m) = 2"^, where m G N, then An is a subset of Z, i.e.. An C Z. In this case, the 
nth particle interacts only with chain particles with numbers n ± 2, n ± 4, n ± 8, n ± 16 .... 

The power law a{m) = 6™, where 6 G N, and 6 > 1, can be realized for compact structure 
of linear polymer molecules. For example, a linear polymer molecule is not a straight line. 
Usually this molecule can be considered as a compact object. It is well-known that "tertiary 
structure" of proteins refers to the overall folding of the entire polypeptide chain into a specific 
3D shape fiQ[ \T7\ [T8] . The tertiary structure of enzymes is often a compact, globular shape 
P!B|[T7]. In this case, we can consider that the chain particle is interacted with particles inside a 
sphere with radius R. Then only some subset of chain particles act on nth particle. We assume 
that nth particle is interacted only with kth particles with k = n ± a{m), where a{m) G N 
and m = 1,2,3,.... The polymer can be a mass fractal object [19j. For fractal compactified 
linear polymer chains, we have the power-law N{R) ~ R^, where 2 < D < 3 and N{R) is the 
number of chain particles in the ball with radius R. Then we suppose that a{m) is exponential 
type function such that a{m) = 6*", where 6 > 1 and 6 G N. This function defines the fractal 
long-range interaction. 

Let us consider equation ([3]) with a fractal long-range interaction. Using the smooth function 
dS]), we obtain the differential- difference equation 

d^u{x,t) = — 2_. b{m)[u{x + a{m)h,t) — 2u{x,t) +u{x — a{m)h,t)]. (8) 

m=— oo 

Using 

exp{—ia{m)hdx} u{x, t) = u{x + a{m)h, t), (9) 

equation ([8]) can be presented as the pseudo-differential equation 

dMx, t) + ^ Y. ^(^) ^i^' {-'-^^^d,) «(^, t) = 0. (10) 

m=—OD 



For a{m) = 1 and b{m) = 6mo, this equation gives equation ([7]) that describes oscillations for 
the case of the nearest-neighbor interaction. 

Let us consider the pseudo-differential operator 

m=—oo 

The function 

^(x, k) = A exp{ikx) (12) 

is an eigenfunction of this operator: 

C'^{x,k) = X{k)'^{x,k). 
Here A(/c) is the eigenvalue of the operator, such that 

X{k) = 2 2 b{m) sin' (^^^k). 

m=—oo 

Using 2sin^(a;/2) = 1 — cos(q;), we obtain 

+ 00 

A(A;) = Y^ b{m)[l-cos{ha{m)k)]. (13) 

m=— oo 
If 

a{m) = a"", b{m) = a^^-'^"", a en, a > 1, (14) 

then equation (fT3!) has the form 

A(A;) = C{hk), 

where C{z) is the cosine Weierstrass-Mandelbrot function 



C{z)= J2 a^''-'^"^[l-cos(a(^-2)-^)]. 

m=—OD 

The box-counting dimension of the graph of this function is D. The operator (TTTj) can be 
called the Weierstrass-Mandelbrot operator. The spectral graph {k,C{hk)) of this operator is 
a fractal set with dimension D. 

Using the operator (TTT]) . equation (TTUI) takes the form 

d^u{x,t) + -—Cu{x,t) = 0. 



It is not hard to prove that the dispersion law for the chain with the long-range interactions is 
described by equations ([3]) and (IT^ has the form 

uj^ + %C{hk) = 0. 

Then the graph {k,uj{k)) is a fractal. Note that the group velocity Vgroup = duj{k)/dk for the 
plane waves cannot be find, since C{z) is the nowhere differentiable function. 

Let us consider a generalization of conditions (1141) in the form 

a{m) = a'", b{m) = 6™, (15) 

where a G N, 6 < 1, and a > 1. For b = a^^^, we have (IT^ . If we use (TT^ . we can obtain the 
pseudo-differential equation 

dM^, t) + ^ sin^ (-|a.) u{x, t) + ^ £ &- sin^ i''-^^^) "(^' ^)- ^^S) 



m=l 



This equation can be presented in the form 



d^u{x,t) + ^ sm'(-^—d^) u{x,t) + M'c\{x,t) = ^^^V' cos(-tha"'d^u{x,t), (17) 



m=l 



where M^ = 46//i^(l — 6). Equation (IT71) describes the oscillations in the case of the fractal 
LRI. The left-hand side of equation (ITTj) in the limit h —>■ gives the Klein-Gordon equation 

c-'dtu{x, t) - dlu{x, t) + M'u{x, t) = 0. 

The right-hand side of equation flT7|) describes a nonlocal part of the interaction. Note that the 
pseudo-differential operator 

+ 00 

A = ^ 6™ cos (-iha"'d^ , (18) 

m=0 

which is used in equation (IT7I) . has the eigenfunctions flT^ such that the eigenvalues is the 
Weierstrass function W{hk/7i), where W{x) is defined by equation ([1]). The box-counting 
dimension of the graph of the Weierstrass function W{x) is (j2]). The operator (TT8|) can be 
called the Weierstrass operator. The spectral graph {k, W{hk/TT)) of this operator is a fractal 
set with dimension 2 -|- ln(6)/ln(a). 



4 Conclusion 

In this paper, we prove that the chains with long-range interaction can demonstrate fractal 
properties. We consider chains with long-range interactions such that each nth particle is 
interacted only with chain particles with the numbers n ± a{m), where m = 1,2,3, .... The 
exponential type functions a{m) = b"^, where 6 > 1 is integer, are used to define a fractal long- 
range interaction. The equations of chain oscillations are characterized by dispersion laws that 
are represented by Weierstrass and Weierstrass-Mandelbrot (fractal) functions. The suggested 
chains with long-range interactions can be considered as a simple model for linear polymers 
that are compact, fractal globular shape. 
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